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ABSTRACT 

We have constructed models of quasi-periodic variability of X-ray emission from ac- 
creting compact objects. Assuming a general scenario of a propagation model of variability, 
with inverse Compton upscatering as the emission mechanism, we have considered a number 
of cases for the periodic modulation: modulation of the plasma heating rate, cooling rate by 
external soft photons, and the amplitude of the reprocessed component. We have computed 
various observational characteristics which can be compared to good quality data. These in- 
clude Fourier-frequency resolved spectra and results of cross-correlation analysis between 
light-curves at different energies. Each model of modulation predicts specific observational 
signatures, which help in identifying the physical processes driving QPO emission in accret- 
ing sources. 



Key words: accretion, accretion disc 
ries: close - X-rays: binaries 



instabilities - radiation mechanisms: thermal - bina- 



1 INTRODUCTION 

X-ray lightcurves from accreting compact objects (neutron star X- 
ray binaries, black hole X-ray binaries and active galactic nuclei) 
are generally non-periodic, showing significant variability in broad 
range of time scales. Consequently, their power spectra are broad, 
extending through a few decades in Fourier frequency (van der 
Klis 1995; Markowitz et al. 2003). Nevertheless, periodic, or quasi- 
periodic components appear often in the lightcurves of the stellar 
systems. Typically, these quasi-periodic oscillations (QPO) contain 
only a small fraction of total variability power, and are hardly, if at 
all, recognizable in lightcurves in the time domain. Usually, they 
are only clearly visible in power density spectra (PDS), where a 
large number of oscillations is co-added. The observed QPO fre- 
quencies range from Keplerian frequency near a surface of a neu- 
tron star (kilo-Hz QPO in some neutron star X-ray binaries) to be- 
low ~ 0.1 Hz (review in Wijnands 2001). 

Despite the fact that a QPO component is almost never the 
dominant light curve component, their significance for our under- 
standing of the accretion process should not be underestimated. In 
accreting black hole systems QPO features seem to appear pref- 
erentially when the state of the source is changing. For example, 
Rutledge et al. (1999) found correlations between spectral state 
transitions and QPO appearance in soft X-ray transient sources in 
Ginga data. Numerous works were devoted to study such corre- 
lations in good quality RXTE data (e.g., Cui et al. 1999; Kalemci 
et al. 2004; Rodriguez et al. 2004). In many neutron star X-ray 
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binaries many long-lived QPO features were observed and studies 
of their time history revealed such interesting phenomena as satura- 
tion of QPO frequency with mass accretion rate (Kaaret et al. 1999) 
or "parallel tracks" phenomenon (van der Klis 2001 and references 
therein). This latter phenomenon is a hysteresis effect similar to 
that observed also on longer time scales in some black hole X-ray 
binaries (Maccarone & Coppi 2003). QPO offer also a method to 
constrain the mass of the central compact object, provided that their 
frequency is identified with a physical frequency of the system. 

A variety of spectral/timing behaviour, including several types 
of QPO, is shown by the black hole binary (BHB) system GRS 
1915+105. This is a peculiar source, most likely due to its high 
mass accretion rate (Sobolewska & Zycki 2004; Done, Wardzinski 
& Gierlinski 2004). Recently, Miller & Homan (2005) noticed 
changes in the intensity of the Fe Kq fluorescent line at 6.4 keV 
during QPO period in this source. This is a potentially important 
clue to the origin of the periodicity, suggesting that variations oc- 
cur in the relative configuration of the cold reprocessing matter and 
the X-ray source. 

All of these demonstrate the importance of understanding the 
origin of the QPO features. Of course, given the variety of QPO 
components it is rather unlikely that a single mechanism could ex- 
plain all of them, despite correlations between different QPO fre- 
quencies (Belloni, Psaltis & van der Klis 2002). 

QPO models considered so far concentrated on identifying the 
frequencies of oscillations (see, e.g., Lee, Abramowicz & Kluzniak 
2004 and references therein, and a review of QPO models in Psaltis 
2001). Rather little attention was paid to the fact that it is the X-ray 
emission, which is modulated. A notable exception is the paper by 
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Giannios & Spruit (2004), where a model of global oscillations of 
the hot inner flow interacting with outer cool disc was constructed. 
Oscillations are driven by feedback loop between heating the cool 
disc by ions from the hot flow and cooling of the hot disc by soft 
photons from the cool disc. The model predicts that r.m.s. variabil- 
ity of the QPO is highest at lower energies, therefore creating a 
rather soft QPO spectrum. 

In this paper we address the problem of QPO origin from 
a complementary point of view, considering how such variable 
hard X-ray emission may be generated. Assuming that the X-ray 
are produced in the thermal (or hybrid thermal-nonthermal; Coppi 
1999) Comptonization process (review in Zdziarski & Gierlinski 
2004), we construct a number of geometrical/phenomenological 
scenarios of oscillations. The main idea here is that - irrespectively 
of the physical mechanism - QPO may be driven by modulations of 
only a few physical parameters which determine the Comptonized 
spectrum: plasma heating rate, cooling rate by soft photons, ampli- 
tude of the reflected component (feedback between the heating and 
cooling), or a temperature of soft seed photons. It is also possible 
that a QPO is a simple geometrical effect of, for example, periodic 
obscuration of the source of emission. For each such case we pre- 
dict various energy and timing characteristics which could be com- 
pared with observations. In particular, we emphasize the Fourier- 
frequency resolved spectroscopy as an important tool (Miyamoto 
et al. 1991; Revnivtsev, Gilfanov & Churazov 1999; Zycki 2002; 
2003, hereafter Z03). 



2 THE MODEL 

The general scenario for time variability considered in this paper is 
a radial propagation model, as described in details in Z03 and Zycki 
(2004). That is, we assume that X-rays are produced by structures 
propagating radially inwards, towards the central compact object 
(Kotov, Churazov & Gilfanov 2001). Dissipation of gravitational 
energy causes a flare of radiation as the structures are approaching 
the center. The X-rays are produced in inverse Compton process on 
thermal electrons, with the soft seed photons provided by optically 
thick cold accretion disc (review in Done 2002). Evolution of en- 
ergy spectrum during a flare of radiation, from softer to harder, pro- 
duces hard X-ray time lags - commonly observed in these sources 
(Poutanen & Fabian 1999; review in Poutanen 2001). The model 
reproduces also the Fourier-frequency resolved spectra (/-spectra; 
Revnivtsev et al. 1999), which get harder with increasing Fourier 
frequency, at least in low/hard state of BHB (Z03). We employ 
the idea of flare avalanches to generate lightcurves with PDS of a 
doubly-broken power law form (Stern & Svensson 1996; Poutanen 
& Fabian 1999). The correlations between the flares (avalanches) 
go some way towards explaining the flux-cr relation of Uttley & 
McHardy (2001). 

We note that the propagation model was formulated for 
low/hard state of BHB, whose spectral/timing properties are known 
best. QPO are mostly observed in softer/brighter states of BHB. 
However, there is certainly a continuity of spectral/timing proper- 
ties when the sources make transitions between the states (see, e.g., 
Zdziarski & Gierlinski 2004; Done & Gierlinski 2004; Kubota & 
Done 2004), so it does seem reasonable to adopt the same general 
scenario for the soft states as well. On the other hand, data anal- 
yses do show differences between some timing properties in hard 
and soft states, for example in Fourier-frequency resolved spectra 
(Revnivtsev et al. 1999 for hard state; Sobolewska & Zycki 2005 
for soft state). 



The flare of radiation is described by a specific time depen- 
dence of plasma heating rate, expressed as compactness parameter, 
I = (Lj D)(aT /m e c 3 ), where L is the luminosity and D is the 
characteristic dimension of the emission region. We adopt the fol- 
lowing descriptions of lh (Z03), 



lh,o(t) = Ar{t)~ 



1 - 



6 
r(t) 



< t < tn 



(1) 



which describes the dissipation of gravitational energy in an accret- 
ing ring. Initial and final radial positions of the propagating struc- 
tures are denoted r; and rf, respectively. i max is the duration of a 
flare, computed assuming motion of the emission region at a frac- 
tion of free-fall velocity (Z03). 

The cooling of the plasma is described by cooling compact- 
ness, l s . In the basic model the main source of soft photons is re- 
processing of hard X-rays, which is described as (Z03), 



i.(t) = NC(t)l h (t), 



(2) 



where C(i) has the meaning of a covering factor of the reprocess- 
ing matter. It is described by 



C[r(t)] 



[r(t)/r tI ] 



for r(t) > r tr 
for r(i) < r tr , 



(3) 



i.e. it equals 1 when the emitting structure is above the cold disc, 
but it gradually decreases to approaching the black hole. The de- 
cline of C(r) below the truncation radius, r tr , has indeed to be 
rather gradual, as discussed in some detail in Z03; specifically the 
exponent is assumed 7 = 2. 

As indicated above, the idea for generating a QPO is then to 
impose an additional periodic modulation of one or more of the 
above time dependencies. For example, if the QPO were a result of 
oscillations of the accretion disc, then irrespectively of physical ori- 
gin of the oscillations, they could affect the flux of soft photons, or 
the temperature of the soft emission, To, or the amplitude of the re- 
flected component, R, or some combinations of the above. This will 
produce a definite signature in, for example, energy dependence of 
QPO strength, or, related quantity, energy spectra computed at the 
QPO frequency. 

The Comptonized spectra are computed with the EQPAIR code 
(Coppi 1999). The code computes thermal and pair balance of a 
spherical plasma cloud, given the heating rate of the plasma and 
cooling rate by soft photons, and other necessary parameters. It 
then computes the spectrum of emerging radiation. In our case, at 
each time step the heating and cooling compactnesses are calcu- 
lated from eq. and with the additional periodic modulation 
as detailed below, and the spectrum of emergent radiation is com- 
puted. The reprocessed component is then added to it, with relative 
amplitude, Q/2n — C(t), where fi/27r is the solid angle of the 
reprocessor from the X-ray source. Narrow Fe Kq fluorescent line 
is added to the Compton reflected continuum. Its energy is fixed at 
6.4 keV, while its equivalent width (EW) depends of the continuum 
slope (Zycki & Czerny 1994; Z03), roughly covering the range of 
1 to 1.3 keV. 

A sequence of spectra is computed with total duration of 256 
sec. with time step of 1/32 sec. It is then analysed using standard 
tools of Fourier analysis. The /-spectra are computed as described 
in Revnivtsev et al. (1999) and Zycki (2002). An implicit assump- 
tion here is that the radius of the flaring region is constant in time, 
that is, the flare luminosity is represented by the compactness pa- 
rameter. One may imagine more complex models where both D 
and L vary, which may lead to luminosity-spectrum relations dif- 
ferent than those presented in this paper. 
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Figure 1. Energy spectra from a test case of modulation of reflection am- 
plitude only fSec l3.ll . Solid curve is the time averaged spectrum, squares 
(green online) is the r.m.s. spectrum of the broad band noise component, 
while circles (red online) is the r.m.s. spectrum of the QPO. The QPO was 
modeled by varying the amplitude of reflection only, and its r.m.s. spectrum 
reproduces the shape of the reflected component. 

3 RESULTS 

Computations are performed for a 10 M© black hole system. Flare 
duration (travel time), r, is generated from a probability distribu- 
tion P(t) oc t -1 between r m i n = 0.005 sec and r max = 2 
sec, so that the PDS roughly matches PDS of Cyg X-l (Z03). Ini- 
tially, n is set to 100 R g but it is set smaller for short flares, so 
that the propagation velocity does not exceed the free-fall velocity. 
The normalization constant in eq. Q is such that the maximum of 
ih(t) = 100. The parameters are chosen so that the time average 
spectra are rather hard: normalization constant in eq.|2| N = 0.4, 
and r tr = 20i? g , which give final spectrum with V ~ 1.8. The 
QPO frequency is assumed 1 Hz. 

3.1 A test case - modulation of reflection amplitude 

Let us consider first a simple case of modulation of amplitude of 
the reprocessed component, Q/2n. That is, we use un-modified 
eqs. 0, J5J and {3J to describe heating and cooling of the plasma, 
but the reflection amplitude is 

Q/2n(t) = A qpo sm 2 (2nf qpo t + qpo ), (4) 

with j4 qpo = 2 and / qpo = 0.5 Hz. We call it a test case because 
it is not a self-consistent situation in our model, since fl/2n should 
equal to the feedback function C(t). It does, nevertheless, demon- 
strate well the meaning of the r.m.s. spectrum. The r.m.s. spectra 
were obtained by first fitting to the PDS in each energy channel a 
model consisting of a doubly broken power law, 

f f°, if /</M 
P(f) =N P xl Nif- 1 , if /,,,! < / < f b>2 (5) 

{ N 2 f~ 2 , if f>f b ,2 

(where N± — fb,i and N2 = fb,ifb,2 ensure continuity of P (/)), 
with a Lorentzian QPO feature. Fitted parameters are the overall 
normalization Np, two break frequencies fb,i, fb,2 and normal- 
ization and width of the Lorentzian. The rms(i5) are then simply 
computed as integrals over / of the continuum P(f) and the QPO. 



The QPO r.m.s. spectrum plotted in Fig. Q reproduces the energy 
spectrum of the process which generated the QPO. 

3.2 Modulation of the heating rate 

Here we consider periodic modulation of the heating rate, 1^. The 
motivation for this case comes from the description of broad band 
X-ray power spectra as a sum of a number of Lorentzians (e.g., 
Miyamoto et al. 1991; Nowak 2000). If taken seriously, it may be 
interpreted that X-rays (not only QPO) are emitted in the form 
of damped or forced oscillators. One might then imagine that a 
QPO is one of those oscillatory modes, of rather higher quality 
factor than oscillations producing the broad band noise. Further- 
more, Maccarone, Coppi & Poutanen (2000) demonstrated that in 
the low/hard state the flares have to be driven by the heating rate 
rather than cooling flux of soft photons. Specifically, we assume 
then 

lh(t - to) = lh,o(t - to) [1 + -4 qP o sin(27r/ qP oi + 4> qpo )] , (6) 

where to is start time of a given flare. The time t is a global time, 
thus the above formula describes really a background modulation 
acting on a collection of flares, whose start time and durations 
are randomly generated from appropriate probability distributions 
(Z03). If the phase, qpo , was constant, the modulation would be a 
perfectly periodic, thus producing a narrow feature in PDS. In or- 
der to produce a quasi-periodic modulation, we change the phase 
abruptly at random moments in time. Such a procedure is motivated 
by results of Morgan, Remillard & Greiner (1997), who found ran- 
dom walk of QPO phase in GRS 1915+105. Parameters A qpo and 
/qpo are constant. 

Of crucial importance is the response of the soft flux, l s , to the 
modulation of lh- In the basic model of Poutanen & Fabian (1999) 
the soft flux comes only from reprocessing of hard X-rays, which 
means that the former must respond coherently to the variations 
of the latter. This assumption was made so that the model repro- 
duces the spectral slope - reflection amplitude correlation observed 
in low/hard state (Zdziarski, Lubinski & Smith 1999). However, 
in soft states the disc emission usually dominates the hard X-ray 
emission, therefore we have to allow for the possibility that they 
are at least partially independent. Consequently, we consider also 
the situation where the soft flux might respond only to the flare of 
radiation, described by lh,o{t), but it does not respond to the quasi- 
periodic component. 

3.2.1 Cooling rate responding to QPO 

Here we assume that the flux of soft photons cooling the plasma is 
related to the heating as 

l B (t) = NC(t) x l h (t), (7) 

that is, the soft flux comes from reprocessing of the entire flux of 
hard X-rays, including the quasi-periodic modulation. We assume 
that any delay in the response of l s (t) is negligible. This is certainly 
correct for QPO of considered frequency, / ^ 10 Hz, since the 
light travel time through 10J?sch is only 10 -3 s. 

The results are presented in Fig. [2] Panel (a) shows the energy 
spectra folded with the period of the QPO, T qpo = 27r// qpo , and 
divided by the time average spectrum. No clear spectral evolution 
during the QPO is seen. This is obviously a consequence of eq.0 
since the spectral slope depends only on the ratio which is 

unaffected by the QPO. Panel (b) shows the frequency-integrated 
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Figure 2. Simulation results for the case of quasi-periodic modulation of the heating rate, with coherent response of the cooling rate (see Sec. 13.2. IV (a) 
shows the energy spectra folded with the period of the QPO, T qpo = 2iz/f qpo , and divided by the time average spectrum. The spectra show no clear spectral 
evolution during the QPO. (b) shows the frequency-integrated rms variability as a function of energy for both the broad band noise component (excluding the 
QPO; solid squares, green online) and the QPO (circles; red online). The QPO rms is basically independent of energy, as a consequence of the lack of spectral 
variability, while the noise rms increases somewhat with energy, (c) - energy spectra: time averaged (solid curve), the variable component (filled squares; 
green online), (d) shows /-resolved spectra, computed in the following ranges of /: 0.04—0.06 Hz, 0.4-0.7 Hz, 0.9-1.3 Hz (QPO; blue online), 1.7-2.5 Hz, 
4.1-6.1 Hz and 1 1-16 Hz. There is a weak overall trend for hardening of the spectra with / and the QPO spectrum is not much different than for other values 
of frequency, (e) shows the EW of the Fe Ka line as a function of / (thick histogram; red online), comparing it to the basic case of no QPO (thin histogram). 
Again, in the considered case there is no additional signal at / qpo . (f) and (g) show results of cross-correlation analysis of light curves at 3 keV and 9 keV. The 
time lags show a minimum at / qpo , which seems to be the only signature of the QPO. The coherence function (g) is equal to 1 in almost entire range of /. (h) 
shows the PDS from light curve at 9 keV. 




Figure 3. Results for the case of quasi-periodic modulation of the heating rate, with no response of the cooling rate (see Sec. 13.2.21 . The QPO spectrum is 
harder here than the /-spectra at neighboring frequencies. It is also harder than the time averaged spectrum. 
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r.m.s. variability as a function of energy for both the broad band 
noise component (excluding the QPO) and the QPO. This is com- 
puted as explained in Sec. 13. II The QPO rms(_E) is basically inde- 
pendent of energy, again as a consequence of the lack of spectral 
variability. The noise rms(_E) increases somewhat with energy, be- 
cause of the spectral evolution (hardening) during each flare. Panel 
(c) shows energy spectra in E Fe representation: time averaged 
spectrum, r.m.s. spectrum of the variable component and the QPO 
r.m.s. spectrum. The two latter spectra are simply the rms(iJ) de- 
pendencies from previous panel multiplied by time averaged count 
rate Fe- Panel (d) show the /-spectra computed in the following 
ranges of /: 0.04-0.06 Hz, 0.4-0.7 Hz, 0.9-1.3 Hz (QPO), 1.7- 
2.5 Hz, 4.1-6.1 Hz and 11-16 Hz. There is a weak overall trend 
for hardening of the spectra with / and the QPO /-spectrum (i.e. 
the spectrum at / around / qpo ) is not much different than for other 
values of frequency 1 . Panel (e) shows the equivalent width (EW) 
of the Fe Ka line as a function of Fourier frequency, comparing it 
to the basic case of no QPO. Again, in the considered case there 
is no additional signal at / qpo . Panels (f) and (g) show results of 
cross-correlation analysis of light curves at 3 keV and 9 keV. By 
construction, the model produces hard X-ray time lags. The time 
lags show what seems to be the only signature of QPO in the con- 
sidered case, namely the lags show a clear minimum at / qpo . The 
coherence function is equal to 1 in almost entire range of /, in par- 
ticular there is no signature of QPO here. Finally, (h) shows the 
PDS from light curve at 9 keV. 

3.2.2 Cooling rate not responding to QPO 

In this case we assume that the soft flux comes from reprocessing 
of the non-periodic flare of hard X-rays only, 

k(t) = NC(t) x i M (t). (8) 

It is not entirely clear how this situation may be real- 
ized geometrically, since it would probably require a specifically 
anisotropic emission. It is nevertheless an interesting one to exam- 
ine because of clear observational signatures. Results, plotted in 
Fig. show an obvious pattern of spectral variability during the 
QPO: the brighter the spectra the harder they are 2 . A consequence 
of this is that the QPO rms increases with energy, and the QPO en- 
ergy spectrum is much harder than the other /-spectra. The harder 
QPO spectrum produces somewhat larger EW of the Ka line at 
/qpo- Similarly to the previous case, the 9 keV vs. 3 keV timelags 
show some sign changing around / qpo , and the coherence function 
remains equal to 1 at all /. 

3.3 Modulation of the soft flux 

Here we assume that the flux of soft photons is additionally modu- 
lated in quasi-periodic manner. Again, we consider two situations: 
modulation of the flux with constant temperature of soft seed pho- 
tons, and modulation of the flux with corresponding change of the 
soft photons temperature. 

1 The QPO /-spectrum is constructed from total P(f) around / qpo , i.e., 
containing contribution both from the broad band noise and the QPO, while 
the previously considered QPO r.m.s. spectrum was constructed from the 
component describing the QPO only 

2 Note that our model produces a broad band PDS with a QPO, which 
means that individual spectra may not obey the trends; the trends are visible 
as the shapes of envelopes of QPO period-folded spectra. 



3.3. 1 Constant soft photons temperature 

We assume that the temperature of the soft seed photons is con- 
stant in time, while their flux is additionally modulated. Specifi- 
cally, £h(t) = hi,o(t) (eq-E) and 

ls(t - to) = NC(t - t Q )l h (t) + Ajp [sin(27r/ qpo f. + <?Vo)] 2 , ( 9 ) 

where the square of the sine function assures that l s > 0. Such 
form of l s could correspond to accretion disc pulsations modulating 
the intrinsic (rather than reprocessed) flux of soft photons. We now 
assume / qpo = 0.5 Hz, so that the QPO is still at 1 Hz in PDS. 

Spectral variability driven by the soft flux produce a charac- 
teristic pivoting behaviour, with the pivot energy intermediate be- 
tween To and kT c . This is evident in our simulations presented in 
Fig.|4| where the pivot energy is about 20 keV. As a consequence, 
the rms(_E) has a broad minimum around 20 keV. The QPO spec- 
trum is in consequence much softer than the other /-spectra, at least 
up to ~ 10 ke V. This in turn explains the small EW of Ka line at the 
QPO frequency. Time lags between 3 keV and 9 keV light curves 
do not show any clear feature at / qpo , but a sharp drop of the coher- 
ence function at / qpo is observed. This is caused by the existence of 
two independent Fourier components at / qpo : one from the broad 
band noise, the other from the QPO. The lack of QPO in PDS from 
the 9 keV light curve is consistent with the presence of minimum 
in rms(£). 



3.3.2 Modulation of the seed photons temperature 

We assume that the temperature of the soft photons, Tsoft , is mod- 
ulated as 

T soft (t) = To [1 + A qpo sin(27r/qp t + </> qpo )] , (10) 
and the corresponding modulation of the soft flux is given by 

l^t) = NC(t)l h (t) + N(^y , (11) 

where we assume A qpo = 0.1 and / qpo = 1 Hz. The above form 
describes two sources of soft photons: in addition to the reprocessed 
photons, there is an additional flux from periodically modulated 
intrinsic dissipation. We assume that the geometrical dependence 
of the modulated flux is different than that of the reprocessed flux 
and that the covering factor C(t) does not appear in the second 
term in the expression for cooling compactness. 

The results (Fig[3J are rather different than those for previ- 
ous case of constant T so f t . Strongest variability is now at low ener- 
gies, around E = fcT ao ft, and its amplitude decreases with energy. 
Consequently, the QPO spectra are soft and the Fe Ka line EW is 
lower than our reference case. It does not show any clear feature at 
/q po . Time lags do show some sign changes around / qpo , while the 
coherence shows a complex behaviour, with local minima around 
/qpo, but a maximum equal to 1 at / qpo . There is also a harmonic 
QPO at 2/qpo, which is related to the fact that the flux modulation 
is not a simple sinusoidal modulation. The harmonic is strong only 
at the energies where the soft flux appears. 

3.4 Modulation of the covering factor of the reprocessor 

Here we assume that the covering factor of cold matter is addition- 
ally quasi-periodically modulated, 

C(t) = Co (£)[! + A qpo sm(27r/q P ot + qP o)], (12) 
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Figure 4. Results for the case of quasi-periodic modulation of the soft cooling flux (see Sec. l3.3.fl with constant temperature of the seed photons. See caption 
to Fig.|2]for detailed description of panels and curves. Characteristic pivoting behaviour resulting from modulations of l s gives broad minimum in rms(i?) at 
~ 10 keV and, consequently, softer QPO spectra and a minimum in Ka line EW(/) at / = / qpo - 



with ^qpo = 0.5, /qpo = 1 Hz, and qpo generated as described 
earlier. Both the heating rate and cooling rate are described by 
their basic equations, eq. \\\ and eq. [2] respectively. Compared to 
the previous case of modulation of l s , now the amplitude of the 
reprocessed component (Compton reflected continuum with con- 
sistent Fe Ka line) varies as well, Q/2ir(t) = C(t). This scenario 
may correspond to disc warping due to, e.g., Lens-Thirring pre- 
cession (e.g., Miller & Homan 2005). Results of computations are 
presented in Fig.|6| The QPO-folded energy spectra show clear piv- 
oting around 30-50 keV. The variability of the Ka line with the 
QPO phase is noticeable even in those folded spectra. The QPO 
rms(iJ) dependence shows an increase of variability rms at the line 
energy, and a broad minimum around 20-50 keV. The QPO spec- 
trum is then softer than the other /-spectra, while the EW of the 
Ka line has a sharp maximum at / qpo - This is a clear difference 
with the previous case of l s modulation, where a soft QPO spec- 
trum was accompanied by a weak line. Here, the correlation be- 
tween the spectral slope and amplitude of reflection is obeyed in 
the /-spectra, because of eq.|2| Thus a soft continuum is accompa- 
nied by a strong reflection component, hence strong Ka line. The 3 
vs. 9 keV time lags appear to reverse sign at / qpo (but are positive 
on both sides of it), but there is only very weak loss of coherence 
between the two light curves at / qpo - 



3.5 Modulation of absorption column 

One might imagine that modulation of the observed flux could be 
produced simply as a result of obscuration by an absorber. Here 
the results fall into either of two extreme categories. If the variable 
obscuration affects the source spectrum, the results are very com- 
plex. For example, the power spectra contain a number of strong 
harmonics at 2x, 3x, ...of the fundamental frequency. Also the 
coherence function then shows a number of minima at the same 



/'s. This is obviously because the way the spectra are distorted by 
absorption is not linear in TVh, thus the complexity appears even if 
Nh is modulated periodically. On the other hand, if the absorber 
is very optically thick tt 3> 5, so that practically no photons are 
transmitted and the only variable factor is the covering factor, then 
only the luminosity is modulated, but no spectral signatures are ob- 
served. 



4 DISCUSSION 

We have constructed models of X-ray QPO in accreting compact 
systems, and computed specific energy-spectral properties of the 
QPO in each model. We have assumed general radial propagation 
model, where the X-ray emission is produced by structures moving 
radially towards the central black hole (Kotov et al. 2001 ; Z03). As- 
suming that the relevant radiative process is thermal Comptoniza- 
tion, we considered cases when the quasi-periodically modulated 
factor is plasma heating rate, cooling rate by soft photons, ampli- 
tude of the reprocessed component (feedback factor between the 
heating rate and cooling rate), or the temperature of soft photons. 
For each case we have constructed a set of observable quantities, 
which may be compared to good quality data. These are generally 
of two types: energy dependencies at various Fourier frequencies, 
and Fourier frequency dependencies, showing possible signatures 
at /qpo- QPO amplitude as a function of E is the obvious charac- 
teristic, but it is really a powerful tool for constraining models only 
when used in conjunction with similar dependencies at other /. 
Such Fourier spectra show distinct signatures at / qpo in most con- 
sidered models. The QPO spectrum is harder than other /-spectra 
when lh is modulated and l s does not respond fully to that modula- 
tion. The QPO spectrum is softer than other /-spectra if the mod- 
ulated factor is l s or the reflection amplitude, Q,/2ir (assuming the 
latter controls the feedback between lh and l s ). The difference be- 
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Figure 5. Results for the case of quasi-periodic modulation of the soft cooling flux, with corresponding change of seed photons temperature (Sec. l3.3l2l . 
are quite different than those with constant T so ft fFig.HI. Maximum variability is now at low energies, and the QPO spectrum is very soft. Importantly, this 
modulation produces also a harmonic of the QPO, although only at low energies (it is no longer visible at 9 keV). Crosses in panel (c) (blue online) show 
the spectrum of the QPO harmonic. Coherence function (g) shows very complex behaviour: it seems to have a maxima at / qpo and 2/ qpo , with values much 
lower than 1 around these two frequencies. 



tween modulation of l s and Q,/2n is in the behaviour of Fe Ka 
line as a function of /: for l s modulation it shows local minimum 
at /qpo, while it shows local maximum at / qpo for Q/2n modula- 
tion. The coherence function also differs between the case of l s and 
Q./2-k modulation: it shows a dip at / qpo (weaker coherence) in the 
former case only. 

The significance of the presented models is in offering possi- 
bilities of providing hints as to the physical models of QPO. The 
basic question is whether the modulation occurs as a result of in- 
stabilities in the hot plasma or in the cold disc. The former would 
produce modulation of the heating rate, while the latter modulation 
of the cooling rate, soft photons temperature, reflection amplitude, 
or some combination of these. Candidates for the cold disc instabil- 
ities have been discussed in literature for some time now (review in 
Kato, Fukue & Mineshige 1998), but no work was done on conse- 
quences of these on X-ray spectra. This is of course not surprising 
because physical models of transitions between and interactions of 
the cold disc and hot corona are at a very early stages of devel- 
opment (Rozariska & Czerny 2000, Giannios & Spruit 2004 and 
references therein). 

Models presented in this paper certainly do not cover all pos- 
sibilities. One can imagine other, more complex scenarios. In par- 
ticular, we have assumed here purely thermal Comptonization, 
while most spectra of sources in bright soft states require hybrid 
thermal/non-thermal plasmas (Gierlinski & Done 2003 and refer- 
ences therein). Fits to such spectra usually imply that both direct 
heating and injection of energetic non-thermal particles contribute 
to the energy input. One can then imagine that a modulation might 
affect one of those channels, or their relative importance. We dis- 
cuss some of such models in Sobolewska & Zycki (2005) as a 
possible description of low frequency QPO in XTE J1550-564 in 
anomalous very high state. 



Generally, observational data on QPO energy dependencies 
from a number of black hole binaries point to a complex picture 
(Sobolewska & Zycki 2005), although certain correlations do ap- 
pear. The QPO r.m.s. spectra of sources in hard states seem to be 
softer than the time averaged spectra, while the QPO spectra of 
sources in soft state appear to be harder than time averaged spectra. 
According to our analysis, this would point to driving the modu- 
lation by plasma heating in soft state, while in the hard state the 
driver would be related to changes in the cold disc. This would cor- 
respond to the modes of long-term variability discussed for Cyg 
X-l (Zdziarski et al. 2002), perhaps suggesting that the low fre- 
quency QPO (which appear preferentially during state transitions; 
Rutledge et al. 1999), are simply faster and more coherent version 
of the general variability pattern operating in a given spectral state. 
We discuss this in more detail in Sobolewska Zycki (2005). 



5 CONCLUSIONS 

Signatures of quasi-periodic modulation of physical parameters in 
spectral and timing observables are complex, depending on details 
of a particular model. The QPO r.m.s. spectra can be harder than 
time average spectra if the the plasma heating rate is modulated. 
The QPO r.m.s. spectra can be softer than the time average spec- 
trum if parameters of the cool disc are modulated: flux of cooling 
soft photons or covering factor of the cold plasma. 
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Figure 6. Results for the case of quasi-periodic modulation of the covering factor of the cold reprocessing matter, C'(t), (see Sec. 13. 41 . which controls both 
the feedback between heating and cooling of the hot plasma (eq.|2j, and the amplitude of reflection, Q/2w. Strong Fe Ka line present in the QPO spectra is 
the most characteristic feature of this model. 
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